Fermi surface truncation from thermal nematic fluctuations 
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We analyze how thermal fluctuations near a finite temperature nematic phase transition affect the 
spectral function j4(k, cj) for single-electron excitations in a two-dimensional metal. Perturbation 
theory yields a splitting of the quasi-particle peak with a d-wave form factor, reminiscent of a 
pseudogap. We present a resummation of contributions to all orders in the Gaussian fluctuation 
regime. Instead of a splitting, the resulting spectral function exhibits a pronounced broadening of the 
quasi-particle peak, which varies strongly around the Fermi surface and vanishes upon approaching 
the Brillouin zone diagonal. The Fermi surface obtained from a Brillouin zone plot of A(k, 0) seems 
truncated to Fermi arcs. 
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The concept of nematic order in interacting electron 
liquids has attracted considerable interest over the last 
decade, mostly in the context of (quasi) two-dimensional 
systems [H, 01 ■ In a nematic state an orientational sym- 
metry of the system is spontaneously broken, without 
breaking however the translation invariance. One route 
toward a nematic state is via partial melting of stripe 
order in a doped antiferromagnetic Mott insulator [3|. 
Alternatively, a nematic state can be obtained from a 
Pomeranchuk [3| instability generated by forward scat- 
tering interactions in a normal metal |g| . On a square 
lattice, the most natural candidate for a Pomeranchuk 
instability has dx2^y2 symmetry. 

Signatures of nematic order with a dx2-y2 symme- 
try have been observed in several strongly interacting 
electron materials. A nematic phase with a sharply de- 
fined phase boundary has been established for Sr3Ru2 07 
in a strong magnetic field Nematic order has also 
been observed in the high temperature superconduc- 
tor YBa2Cu30y in transport experiments [8| and neu- 
tron scattering Due to the slight orthorhombic- 
ity of the Cu02 planes one cannot expect a sharp ne- 
matic phase transition in YBa2Cu30y. However, the 
strong temperature dependence of the observed in-plane 
anisotropy indicates that the system develops an intrin- 
sic electronic nematicity, which drastically enhances the 
in-plane anisotropy imposed by the structure l3, 

Nematic fiuctuations close to a continuous nematic 
quantum phase transition naturally lead to non-Fermi 
liquid behavior [l2-14|. For a d-wave Pomeranchuk in- 
stability on a square lattice, the decay rate of electronic 
excitations is strongly momentum dependent along the 
Fermi surface. At the quantum critical point, the de- 
cay rate for single-particle excitations is proportional to 



dj^lwp/"^, where ui is the excitation energy and dk is a form 



factor with d-wave symmetry, such as dk = cos kx — cos ky 



I3I Landau quasi-particles are thus unstable every- 
where on the Fermi surface except at the Brillouin zone 
diagonals, where dk vanishes. The temperature depen- 



dence of the decay rate in the quantum critical regime 
near the quantum critical poin t also differs strongly from 
Fermi liquid behavior [I6, 17 1. 

Previous works on non-Fermi liquid behavior caused 
by nematic fluctuations focussed on the quantum critical 
point and the quantum critical regime at finite temper- 
ature. In this paper we analyze the spectral function 
for single-particle excitations in the thermal fluctuation 
regime near a nematic phase transition at finite temper- 
atures. We show that a perturbative calculation of the 
self-energy in that regime leads to a splitting of the quasi- 
particle peak with a d-wave form factor in the single- 
particle excitation spectrum, reminiscent of a pseudogap. 
However, in a self-consistent calculation the split peak in 
the spectral function is replaced by a single broad peak. 
In the Gaussian fluctuation regime, a summation of ver- 
tex corrections to all orders is possible and confirms the 
self-consistent result. The Fermi surface obtained from 
the peak of the spectral function at zero excitation en- 
ergy is thus smeared by a smooth broadening, which is 
most pronounced near the points (tt, 0) and (0, n) of the 
Brillouin zone, while it gradually decreases toward the 
Brillouin zone diagonal. 

We consider a one-band system of electrons on a square 
lattice with a tight-binding dispersion ek and an effective 
interaction of the form 11311 



1 



Hi = ^^5(q)"d(q) nd{- 
q 



(1) 



where ^^(q) = J2k,a Ck_q/2,o-'^k+q/2,<T are d-wave den- 
sity fiuctuation operators, and L is the number of sites. 
The function g{q) is negative and peaked at q = 0, so 
that forward scattering dominates. An effective interac- 
tion of the form Hi can be obtained from microscopic 
models such as the Hubbard or t-J model 

For sufficiently negative values oi g = g{Q) the interac- 
tion generates a d-wave Pomeranchuk instability leading 
to a nematic state with a spontaneously broken orienta- 
tion symmetry 13, 18-23|. A suitable order parameter 



2 



characterizing the symmetry breaking is provided by the 
expectation value {nd{Q)). Close to the transition (if con- 
tinuous), strong d-wave density fluctuations with a long 
wavelength develop, which lead to a singular effective in- 
teraction. In the quantum critical regime the effective 
interaction is dynamical and of the form [l3, 16 1 



£'kk'(q, I'n) = 



g dkrfk' 



(WeP + eo'|qP + l^«|/Hq|) ' 



(2) 



where — 2'imT is a bosonic Matsubara frequency; ^ is 
the nematic correlation length, while and u are non- 
universal parameters determined by the momentum de- 
pendence of 5(q) and the band structure. In the ther- 
mal fluctuation regime near the finite temperature phase 
transition, quantum [vn ^ 0) fluctuations are cut off by 
temperature, such that only the classical part of the ef- 
fective interaction. 



-Dkk'(q) = £'kk'(q,0) = 



(3) 



with g — g/^Q, is important. 

The nematic transition on a square lattice belongs to 
the two-dimensional Ising universality class. A thermal 
phase transition at a critical temperature Tc > is pos- 
sible, since the dimensionality of the system is above the 
lower critical dimension (one). The correlation length ^ 
diverges at Tc. Approaching the critical temperature, one 
first passes through a Gaussian fluctuation regime, where 
order parameter interactions are not important. Moving 
closer to Tc, one enters the Ginzburg region, where order 
parameter interactions become relevant, and the fiuctua- 
tion propagator acquires an anomalous scaling dimension 



21[. Close to the quantum critical point, the width of 
the Ginzburg region is of order Tc/\ logTd (2^ . 

The momentum resolved spectral function for single- 
particle excitations can be written as 

AiKu.) = -ilmG(k,a.) = -^I-, _ _ _ s(k,a.) 

(4) 

where G'(k, oj) and S](k, w) are the retarded Green func- 
tion and self-energy, respectively. We first compute the 
self-energy perturbatively to first order in the effective 
interaction. The contribution from classical fluctuations 
is given by [l^ 



S(k, ^) = -tJ i?kk(q) G(k - q, c.) 



(5) 



In a non-selfconsistent evaluation of Eq. ([5]) one ap- 
proximates G by the non-interacting Green function 
Go(k, w) — [uj — (ck — m) + iO~^]~^. The self-energy can 
then be computed analytically. The imaginary part has 
been obtained already previously flQj . For momenta close 
to the Fermi surface and small frequencies one finds 



ImE(k, w) 



4wk 



(6) 



where Vk = |Vek| is the velocity of the electrons, k = 
[uj — (ek — fJ')]^/vk, and 1{k) = (1 + k^)~^/^. We assume 
that the Fermi surface does not cross van Hove points, 
such that Wk is finite. The real part of the self energy is 
obtained either by a direct evaluation of Eq. ([SJ or by a 
Kramers-Kronig transformation of the imaginary part as 



ReS(k,w) 



1 — kI(k) 



1 + kI{k) 



(7) 



In Fig. 1 we show results for the spectral function as 
obtained from the non-selfconsistent first order calcula- 
tion of the self-energy. Here and in all further numerical 
results we have chosen a dispersion ek = —2t{coskx + 
cosfcj,) — At' cos kx cos ky — 2t"{cos2kx + cos2ky) with 
hopping amplitudes t = 1, t' = —0.3, and t" = 0.15. 
The lattice constant is one, and the chemical potential 
fi has been chosen such that the electron density is fixed 
aX n — 0.9. The corresponding Fermi surface is closed 
around (tt, tt). The coupling constant \sg = —1.2 and the 
temperature T — 0.15. We have not attempted to com- 
pute the correlation length, since it depends on model 
details such as the momentum dependence of ^(q). In- 
stead we show results for various choices of ^. One can 
see that a pronounced splitting of the quasi-particle peak 
develops for increasing ^, which could be interpreted as 
a fluctuation precursor of the symmetry-broken state. 
However, we now show that the splitting disappears in a 
self-consistent calculation, and it is not restored by vertex 
corrections, at least in the Gaussian fluctuation regime. 




FIG. 1: (Color online) Spectral function j4(k,a;) as obtained 
from a non-selfconsistent perturbative calculation of the self- 
energy for various choices of the correlation length ^. The 
graphs show A{]i,ui) as a function of cj for a fixed momen- 
tum k = (0.69, 2.34) on the Fermi surface remote from the 
Brillouin zone diagonal (see inset). 



In a self-consistent evaluation of Eq. ([5]) , with the inter- 
acting Green function on the right hand side, one has to 
solve an integral equation. The problem can be simplifled 
by decomposing the momentum transfer q in components 
qr and qt normal and tangential to the Fermi surface, re- 
spectively 16[ . The dependence on qt can be neglected in 
the momentum argument of G, such that the gj-integral 
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acts only on the fluctuation propagator Dkk(q), yielding 



-Dk(<7r) 



dqt 
2tt 



i?kk(q) 



~9dl 



(8) 



The momentum dependence of the self-energy I](k, oj) 
can be parametrized by the Fermi momentum k^? closest 
to k and the oriented distance from the Fermi surface 
kr, which carries the sign of ^k- One then obtains the 
one-dimensional integral equation 



to) 



gdi^T f dq, 



An 



Vkpikr — qr) — 1^ + S(fcr — 9r,'^) 



(9) 



Note that S does not depend on kr and w independently, 
but only on the difference w — Vk^kr. The dependence of 
S on ki? enters only parametrically via Wk^ ^^'^ and 
has not been written explicitly. The integral equation © 
can be solved numerically. The results for A(k, ui) differ 
strongly from those suggested by non-selfconsistent per- 
turbation theory. The quasi-particle splitting observed in 
the perturbative calculation (Fig. 1) is wiped out com- 
pletely by the self-energy feedback into G. The spectral 
function exhibits only a single peak with a maximum at 
uj — v\^pkr, even for a very large correlation length ^. 

We now consider higher order contributions not con- 
tained in the self-consistent one- loop approximation ([5]). 
The sum over all self-energy contributions generated by 
thermal fluctuations can be written as 



S(k,L.) 



-T 



<fq 



(27r)2 

X A(k-q/2,w;q,0) 



-Dkk(q)G(k-q,cj) 



(10) 



where A(k, oj; q, v) is the irreducible charge vertex includ- 
ing all vertex corrections. We exploit the fact that domi- 
nant contributions are due to small momentum transfers 
q of order to resum vertex corrections via an asymp- 



totic Ward identity [23j, |2J] . For small q, the charge ver- 
tex is related to the current vertex A and the propagator 
via the Ward identity z/A(k, a;; q, ;/) — q • A(k, w; q, v) = 
G-i(k + q/2,cj -I- v/2) - G-'^{k - q/2,a; - i^/2). In the 
Gaussian fluctuation regime, interactions between order 
parameter fluctuations are not important. In a diagram- 
matic representation of perturbation theory, these in- 
teractions are generated by fermionic loops with more 
than two vertices. Neglecting Feynman diagrams with 
such loops leads to two simplifications. First, the ef- 
fective interaction ([3]) remains unrenormalized. Second, 
diagrams contributing to the charge and current vertices 
involve only an open fermionic line. Since contributions 
with small q dominate, the electron velocity Vk enter- 
ing the current operator is almost conserved such that 
the current vertex can be expressed by the charge ver- 
tex as A (k, oj; q, ;/) = VkA(k, w; q, z^). The latter relation 



holds for each Feynman diagram without fermionic loops. 
Combining this with the Ward identity one obtains, in 
the static limit 1^ = 0, 

G-i(k - q/2, uj) - G-\k + q/2, w) 



Vk • q 



A(k,w;q,0) = 

Vk • u 

(11) 

Inserting Eq. (|lip into Eq. pUj) . and using the Dyson 
equation G~^ = Gq ^ — E, one obtains a closed system of 
equations for E and G. Decomposing the momenta k and 
q in radial and tangential components, and integrating 
£'kk(q) over qt as before (self-consistent solution), one 
finds a one-dimensional linear integral equation for G, 



Vkp kr - 

1 + T 



iO+)G{kr,Lo) 

dqr Dkpjqr) 
2ti vi^^qr 



G{kr-qr,i0) , (12) 



with Di^p{qr) from Eq. dH). 

Note that vertex corrections cannot be summed by the 
above method at a nematic quantum critical point or for 
the related problem of non-relativistic fermions coupled 
to a U{1) gauge field. This is because in these cases 
the dominant momentum transfers are almost tangential 
to the Fermi surface, so that the term Vk • q becomes 
subleading compared to contributions originating from 
fiuctuations of the electron velocity 24 1 . 

The integral equation ([T2|) can be converted to a linear 
differential equation by a Fourier transformation. The 
differential equation can be solved by standard methods. 
The result for the spectral function reads 



A(k,^) 



dxA{x) e'('^-^^F''r-)x/v^, 



(13) 



where 



Aix) 



1 



2TTVl^f 



exp 



dx' 



dx"T 



~9dlp 
2nvl 



■Koix"/0 



/o Jo 

(14) 

Kq is a modified Bessel function. Note that A(k, w) de- 
pends on kr and w only via the difference uj — Uk^^r- 

In Fig. 2 we plot the spectral function A(k, oj) for the 
same parameters as in Fig. 1. The function exhibits only 
a single peak with no trace of a splitting. The splitting 
present in Fig. 1 is therefore an artefact of the perturba- 
tion expansion, at least in the Gaussian regime. Vertex 
corrections do not change the self-consistent one-loop re- 
sult qualitatively. Quantitatively they tend to sharpen 
the peak in A(k, w), but only moderately. For a large 
correlation length ^ the width of the peak in A(k, ui) is 
proportional to -^/log^, corresponding to a peak in the 
imaginary part of the self-energy ImE(k,Lij) c>c ^/k>g^ at 
UJ = Vkpkr- The width of the peak in ImE is also pro- 
portional to y/\og^ and therefore increases with ^. This 
is very different from the perturbative result for ImE, 
Eq. ([B]), where the height of the peak increases rapidly 
with ^, while its width shrinks. 



FIG. 2: (Color online) Spectral function A(k,ijj) as obtained 
from a non-perturbative resummation of contributions from 
thermal fluctuations, including vertex corrections. The choice 
of k and the model parameters are the same as in Fig. 1. 



FIG. 3: (Color online) Intensity plot of the spectral function 
yl(k, 0) as a function of k in the first quarter of the Brillouin 
zone at T = 0.2, for a correlation length ^ — 15. The other 
parameters are the same as in Figs. 1 and 2. 



In the quantum critical regime studied previously 16 1 
the spectral function also exhibits a single peak with a 
temperature dependent broadening. In that regime the 
width of the peak is proportional to T^, with a correlation 
length ^ diverging as (Tj logTj)^^/^ upon approaching 
the quantum critical point at T = 0. 

It is striking that perturbation theory indicates a fluc- 
tuation precursor of the symmetry broken state at lead- 
ing order, which is however not robust with respect to 
higher order contributions. One may compare with the 
case of a charge density wave with a finite wave vec- 
tor, where symmetry breaking opens a gap. Perturbation 
theory indicates a pseudogap above the transition tem- 
perature in such systems, for example for the flux order 
studied in Ref. [25|. It would be interesting to analyze 
the fate of the pseudogap in such systems in a calculation 
beyond perturbation theory. 

Due to the d-wave form factor in the effective inter- 
action, the broadening of the spectral function varies 
strongly in momentum space. To illustrate this, we plot 
an intensity map of A(k, 0) in the first quarter of the Bril- 
louin zone in Fig. 3. A(k, 0) as obtained from Eq. (fT3|) 
diverges on the Brillouin zone diagonal, since dk vanishes 
there. In a more complete model, where other interac- 
tion channels should be added, this divergence will be 
cut off at least by regular (Fermi liquid) contributions 
to ImS of order T^. We have therefore included such a 
regular contribution (— T^). Due to the rapid increase of 
the broadening of the peak in ^(k, 0) upon moving away 
from the Brillouin zone diagonal, and the correspond- 
ing decrease in the peak height, the Fermi surface seems 
truncated to arcs, albeit with fuzzy ends. 

In summary, we have computed the spectral function 
A(k, w) for single electron excitations in the presence of 
critical fluctuations near a thermal nematic phase tran- 
sition in a two-dimensional metal. Leading order pertur- 
bation theory indicates a splitting of the quasi-particle 
peak. However, a resummation of contributions to all 
orders reveals that the splitting is an artefact of pertur- 



bation theory, at least in the Gaussian fluctuation regime. 
The spectral function exhibits a pronounced broadening 
with a d-wave form factor, leading to features reminis- 
cent of Fermi arcs in the Brillouin zone. The qualitative 
shape of the spectral function does not depend on the 
specific choice of parameters. Away from the Brillouin 
zone diagonal, the imaginary part of the self-energy has 
a peak at w = ek — M, in clear contrast to the conventional 
Fermi liquid form. 

It is remarkable that the effect of Gaussian thermal 
fluctuations on electronic excitations could be treated 
non-perturbatively. The method used to sum contribu- 
tions from thermal fluctuations to all orders is not re- 
stricted to nematic fluctuations, but could be applied 
equally well to other critical thermal fluctuations with 
a small wave vector, for example, close to a structural 
phase transition. 

A continuous finite temperature phase transition is 
well established at the roof of the nematic dome found 
for Sr3Ru207 in a strong magnetic field Q- That system 
thus provides an opportunity to observe the fluctuation 
effects computed in this work, but also an experimental 
challenge, since the standard tool to measure the momen- 
tum resolved spectral function, that is, photoemission, is 
hampered by the magnetic field. 

Fermi arcs have been observed in photoemission mea- 
surements of the spectral function A(k.,uj) in various 
high-Tc cuprate compounds 26]. However, the Fermi 
surface truncation in these materials is associated with 
a pseudogap formation, while we obtain only a strongly 
momentum dependent broadening of the spectral func- 
tion. Although experiments indicate electronic nematic- 
ity at least in some cuprates [ij, our results thus show 
that another mechanism needs to be invoked to explain 
the photoemission data. 
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